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Abstract 

By using a path information denned for the measure of the uncer- 
tainty of instable dynamics, a theoretical derivation of Fourier's law of 
heat conduction is given on the basis of maximum information method 
associated with the principle of least action. 

PACS numbers : 05. 45. -a (Nonlinear dynamics); 66.10.Cb (Diffusion); 
05.60. Cd (Classical transport); 05.40.Jc (Brownian motion) 

1 Introduction 

The Fourier's law of heat conduction, very well tested in experiments with 
solids, liquids and gases, is given by 

J(x) = -kV x T{x) (1) 

where J(x) is the heat flux, T(x) the temperature at a position x (in general 
a vector) of the ordinary space and k the thermal conductivity. This law 
is valid not only for the case of steady process, but also for the case where 
temperature varies in time. The theoretical derivation of this experimental 
law for fluids and solids is a long history and thus far remains a challenge for 
theorists of nonequilibrium thermostatisticsfl]. Most of the past efforts for 
deriving this law were focalized on special models like harmonic or anhar- 
monic crystals at steady state and local equilibrium[ ]. A general derivation 
is still missing. 

In this work, we will try to give a generic derivation of this law based on 
maximum information principle. The information we address in this work 
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is a quantification of the uncertainty of dynamical process. It is given by 
Shannon formula [3] 

H=-jr Pk \np k (2) 

fc=i 

with respect to certain probability distribution p k of that process and the 
index k is summed over all the possible evolutions of the system of interest 
in its phase space T. A phase volume Q occupied by the system in T space 
can be partitioned into v cells of volume Si with i = 1,2, ...v in such a way 
that Si fl Sj = (i 7^ j) and U" =1 Sj = Q. A state of the system can be 
represented by a sufficiently small phase cell in coarse graining way. The 
movement of a dynamical system is represented by its trajectories (in the 
sense of classical mechanics) in T space. 

When we look at a nonequilibrium system leaving an initial cell a in the 
T-space for some destinations, we can say that, if the motion of the system 
is regular, there will be only one possible trajectory from a to a final cell 
b, or in other words, there will be only a fine bundle of paths which track 
each other between the initial and the final cells. These trajectories must 
minimize action according to the principle of least action and have unitary 
probability. Any other path should have zero probability. But if the dynamics 
is instable with strong sensitivity to initial condition, the things are different. 
Two points indistinguishable in the initial cell can separate from each other 
exponentially. So from a given initial cell, there may be many possible final 
cells each having a probability to be visited by the system. More than one 
paths between two cells are possible (see Figure 1). 

In this work, the above uncertainty is represented by a path information. 
This path information will be proved to take its maximum when the most 
probable paths are just the paths of least action. Then with the help of the 
model of Brownian motion, the Fick's laws for diffusion and the Fourier law 
for heat conduction are derived. 

2 Maximum path information 

Let us consider an ensemble of large N identical systems leaving the initial 
cell a for some destinations in the phase volume B formed by the final phase 
points occupied by the systems. The travelling time is t ao = h — t a . After 
t a b, all the phase points occupied by the systems are found in the volume B 
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partitioned into cells labelled by b. We observe N kb systems travelling along 
a path kb leading to certain cell b. A path probability can be defined by 
Pk b \a = Nk b /N which is normalized by 

overB w b 

Y Y PhW = 1- (3) 

b k b =l 

where Wb is the number of possible paths from a to a given cell b of the 
volume B. We always suppose each path is characterized by its action A kb \ a 
defined for classical mechanical systems by 

A kb \ a = f L kb (t)dt (4) 



ab 



where L kb (t) is the Lagrangian of the system at time t along the path kb- 
The average action is given by 

overB w b 

A " = Y Y Pk b \aA kbla (5) 
b k b =l 

The uncertainty concerning the choice of paths and final cells by the systems 
is measured by the following Shannon information 

overB 

H a = ~ Y YPkb\a^Pk b \a (6) 
b k b 

which we shall maximized under the constraints associated with Eq.(3) and 
Eq.(5) as follows 



overB w b overB w b 

SiHa + aY, YPk b \a + vY Y Pk b \aA kb \ a ) = 0. (7) 
b k b =l b k b =l 



This leads to 



Pk b \a = ^ exp[-r]A kbla ], (8) 



where the partition function 

overB 



Z a= Y Y e M-VA kb \a]- (9) 



In path integral language, Z a can be given by [4] 

Z a = dpbdxi, / dpdx exp{— 77 / [px — H(p, x)]dt}, (10) 

J overB J a J t a i 

where p = mx is the momentum and H(p, x) is the hamiltonian of the system. 
Here the Lagrangian is given by L kb (t) = px — H(p, x). 
Putting Eq.(8) back into Eq.(6), we get 

d 

H a = \nZ a + r\A a = lnZ a - n— \nZ a . (11) 

It is proved that [5] the distribution Eq.(8) is stable with respect to the 
fluctuation of action. It is also proved that if 77 is positive, Eq.(8) is a least 
action distribution, i.e., the most probable paths are just the paths of least 
action. On the contrary, if 77 is negative, then Eq.(8) is a most action distribu- 
tion which means that the most probable paths maximize action. In any case, 
whatever the sign of the parameter rj, the most probable paths maximizing 
path information always correspond to extremum of action (SA kb \ a = 0). In 
other words, for instable dynamical process, the method of maximum infor- 
mation must be used in order to derive correct probability distributions just as 
the action principle must be used to derive the correct trajectories for regular 
dynamics. 



3 Transition probability of Brownian parti- 
cles 

Now we consider a Markov diffusion process of an ensemble of identical par- 
ticles of mass m idealized by Brownian motion. Suppose a certain path in 
Figure 1 along which a Brownian particle of mass m moves from a to b via a 
path k which is simplified by an intermediate cell k^. Between the three cells 
a, kb, and b situated in real space at x a , Xk b and Xb, respectively, the particle 
is free. The action Ak b \ a of the particle from a to b can be calculated to be 

_ m(x b - x kb ) 2 m(x kb - Xg) 2 
Akbla ~ 2(t b -t kb ) + 2(t ki -Q ' {lZ) 

Then from Eq.(8), we have 

1 

Pk b \a = ^exp 



-mi] 



2{t b 



L k b ) 



^exp 



-mi] 



{x kb -x a y 

2(t kb — t a ) 



(13) 
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The separate normalization of the two factors of this distribution over all 
possible positions x b and all possible paths between a and b gives 



Z kb Z a 



~2{t b -t k y 


d/2 


2(tfc — t a ) 


d/2 


mrj 




mrj 





(14) 



where d is the dimension of the diffusion space. 

On the other hand, according to the solution of diffusion equation[6], the 
transition probability for the particle to go from a to b via kb is 



Pk b \a 



[4nD(t b -t kb )] d /i 
1 

[AnD{t kb -t a )Y^ 



exp 
exp 



4D(t b - 


~t kb ) 


( x k b ~ 


x a ) 2 ' 




-*«). 



(15) 



where D is the diffusion coefficient. A comparison of Eq.(15) with Eq.(13) 
leads to 
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V 



2mD 



(16) 



which implies rj is positive because D > 0. So the highest information distri- 
bution Eq.(8) can be also called least action distribution, i.e., the most prob- 
able paths are just the paths of least action. The inverse of this statement 
is: if the most probable paths minimize action, then the diffusion constant 
D must be positive. 

The physical meaning of r] can be revealed if we use a general relationship 

i 2 



D 



2dr 



6 which leads to 



V 



dr 
ml 2 



(17) 



where / is the mean free path and r the mean free time of the Brownian 
particles. 

4 Fick's laws of diffusion 

With some mathematics, it can be proved that p kb \ a satisfies the following 
dynamical equation: 



—Qf~ = DA xPk b \a 



(18) 
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where is the Laplacian of the diffusion space, x = x b and t = t b . Let 
n a and n b be the particle density at a and b, respectively. The following 
relationship holds generally 

over A 

H b= J2 ^2 U aPk b \a (19) 

a k b 



which is valid for any n a . This means 

dn b 



dt 



DA x n b . (20) 



dn(x,t) 



This is the second Fick's law of diffusion [6]. The first Fick's law J(x) 
—W x nb(x) can be easily derived if we consider matter conservation 
V x ■ J(x,t) where 3(x,t) is the flux of the particle flow. 

The diffusion constant D can be related to the partition function Z a by 
combining Eq.(16) and Eq.(f 1), i.e., 

D _ 1 d_Ag _ 1 d\\nZ a ) 

2mdH a 2m dH a dr] ' 1 ' 



5 Fourier's law of heat conduction 

We consider a crystal idealized by a lattice of identical harmonic oscillators 
each having an energy = N y (x,t)hv where h is the Planck constant, v 
is the frequency of a mode and N u (x,t) is the number of phonons of that 
mode situated at x at time t in x — > x + dx and v — > v + dv. Suppose 
that there is no mass flow and other mode of energy transport in the crystal. 
Heat is transported only through the phonon flow. The phonons of frequency 
v diffuse in the crystal lattice, among the lattice imperfections, impurities 
and other phonons with in addition anharmonic effects[f], just like Brownian 
particles of mass m = hu/c 2 having transition probability Pk b \ a - Let n v = 
N y /dx be the density of phonons which must satisfy 

over A 

n u (x b ,t) = ^Zn u (x a ,t)p kb \ a (22) 

a k b 

and also Eq.(20). 
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The total energy density p(x, t) of phonons at x and time t is given by 

p(x,t) = / hvn u (x,t)dxQ(i>)di>, (23) 
Jo 

where g{v)dvdx is the mode number x — > x + dx and v — > z/ + dv, and z/ m is 
the maximal frequency of the lattice vibration. This implies 

^^ = DA x p(x,t) (24) 

A variation of energy density 5p(x, t) can be related to temperature 
change ST(x,t) by 

Sp{x,t) = c5T(x,t) (25) 

where c is the heat capacity per unit volume supposed constant everywhere 
in the crystal. This leads to 

^ = KA x T(x,t) (26) 

where k = Dc is the heat conductivity. This equation can be recast into 

dT(x,t) 



Of 



KA x T(x,t) (27) 



which should be solved to give the evolution of temperature distribution due 
to the heat flow. When a stationary state is reached, temperature is every- 
where constant, i.e., dT ^ t ^ = A x T(x, t) = 0. So the temperature distribution 
is given by V x T(z) = constant. 

Now considering the energy conservation in an elementary volume be- 
tween x and x + dx in which we have — 8p ^' t ' > = V x ■ J(x, t), Fourier's law of 
heat conduction follows 

J( x ) = -kV x T(x) (28) 

where the time variable t is removed since this law is independent of whether 
J and T vary in time. 
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6 Concluding remarks 



A path information is defined in connection with different possible paths of 
dynamical system moving in its phase space from the initial cell towards 
different final cells. On the basis of the assumption that the paths are char- 
acterized by their actions, we show that the maximum path information leads 
to an exponential path probability distribution of action called least action 
distribution meaning that the most probable paths are just the paths of least 
action. With the help of the model of Brownian motion, we show that, from 
this least action distribution, Fick's laws of diffusion and Fourier's law of 
heat conduction can be derived in a general way without assumptions about 
the state of the transport. 
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Momentum p Final volume B, at time t b 




► 

Position x 



Figure 1: The possible phase space paths for a system to go from a given 
cell a to different cells b of a partition of the final phase volume B during 
the time t a b, each having a probability Pk b \ a - More than one paths from a to 
a given cell b are allowed so we may have = 1, 2...Wb where Wb is the total 
number of paths from a to a cell b. 



9 



